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Abstract 

For stereoscopic systems of imaging atmospheric Cherenkov telescopes (lACTs), a 
key parameter to optimize the sensitivity for VHE 7-ray point sources is the in- 
tertelescope spacing. Using pairs of telescopes of the HEGRA lACT system, the 
sensitivity of two-telescope stereo lACT systems is studied as a function of the tele- 
scope spacing, ranging from 70 m to 140 m. Data taken during the 1997 outburst of 
Mrk 501 are used to evaluate both the detection rates before cuts, and the sensitivity 
for weak signals after cuts to optimize the significance of signals. Detection rates de- 
crease by about 1/3 between the minimum and maximum spacing. The significance 
of signals is essentially independent of distance, in the range investigated. 



1 Introduction 



lACT stereoscopy - the simultaneous observation of air showers with multiple imaging at- 
mospheric Cherenkov telescopes (lACTs) under different viewing angles - has become the 
technique of choice for most of the next-generation instruments for earth-bound gamma-ray 
astronomy in the VHE energy range, such as VERITAS [1] or HESS [2]. The stereoscopic ob- 
servation of air showers allows improved determination of the direction of the primary and 
of its energy, as well as a better suppression of backgrounds, compared to single lACTs. A 
crucial question in the layout of stereoscopic systems of lACTs is the spacing of the identical 
telescopes. Obviously, the spacing of telescopes should be such that at least two telescopes fit 
into the Cherenkov light pool with its typical diameter of around 250 m. For higher energy 
showers the light pool can increase considerably, but at large distances (> 130 m) one observes 
primarily light from shower particles scattered at larger angles. The coincidence rate between 
two telescopes will decrease with increasing spacing; on the other hand, the angle between the 
views and hence the quality of the stereoscopic reconstruction of the shower geometry will im- 
prove. The opinions within the lACT community concerning the optimum geometry of lACT 
arrays differ; the HESS project, for example, initially aimed for a spacing of 100 m between 
adjacent telescopes; more recently, larger distances around 120 m or more are favored [3,4]. 
For the otherwise similar VERITAS array, intertelescope distances of initially 50 m, and later 
80 m to 85 m were foreseen. Among existing systems, the HEGRA 5-telescope lACT system 
has a characteristic spacing of 70 m between its central telescope and the corner telescopes. The 
WHIPPLE-GRANITE two-telescope system [5] had a spacing of 140 m. The Telescope Array 
[6] uses a spacing of 120 m between the first three telescopes, and later 70 m for the full array. 

Optimization of lACT system geometry is heavily based on Monte-Carlo simulations. Over 
the last years, the quality of these simulations has improved significantly, both concerning 
the reliability of the shower simulation [7] and concerning the details of the simulation of the 
telescopes and their readout (see, e.g, [8]). Simulations have been tested extensively, and key 
characteristics, such as the radial distribution of Cherenkov light in the light pool, have been 
verified experimentally [9]. Nevertheless, a direct experimental verification of the dependence 



of the performance of lACT systems on the intertelescope distance would be highly desirable. 

The HEGRA lACT system at the Observatorio del Roque de los Muchachos on La Palma 
consists of five telescopes, four of them arranged roughly in the form of a square with 100 m 
side length, with the fifth telescope in the center. Selecting pairs of telescopes, the performance 
of two-telescope stereo systems can be studied for distances between 70 m (from the central 
telescope to the corner telescopes) and 140 m (across the diagonal of the square), covering 
essentially the entire range of interest. This paper reports the results of such a study, based on 
the large sample of gamma-rays [12] acquired during the 1997 outburst of Mrk 501. 



2 Telescope hardware and data set 



The HEGRA lACT system consists of five telescopes, each with a tessellated mirror of 8.5 m^ 
area and 5 m focal length, and a 271-pixel camera with a 4.3° field of view. The trigger condition 
requires a coincidence of two neighboring pixels above a threshold qo to trigger the camera, and 
a coincidence of triggers from two cameras to record the data. Details of the trigger system 
are given in [10]. Events are reconstructed by parameterizing the images using the Hillas image 
parameters, and by geometrically determining the direction and the impact point of the shower. 
A simple and relatively efficient method for cosmic-ray rejection is based on the difference in the 
width of 7-ray and cosmic-ray images, respectively. Width values are normalized to the mean 
width expected for a gamma-ray image of a given intensity and impact distance, and a mean 
scaled width is calculated by averaging over telescopes. The 7-ray showers exhibit, by definition, 
a mean scaled width around 1, whereas the broader cosmic-ray showers show larger values. A 
cut requiring a m,ean scaled width below 1.2 or 1.3 keeps virtually all gamma-rays and rejects 
a significant fraction of cosmic rays; a cut at 1.0 to 1.1 has lower gamma-ray acceptance, but 
optimizes the significance for the detection of 7-ray sources. The data analysis and performance 
of the system are described in more detail in [11,12]. 

During 1997, when the data for this study were taken, only four of the five telescopes were 
included in the HEGRA lACT system; the fifth telescope (one of the corner telescopes) was 
still equipped with an older camera and was operated in stand-alone mode. The four telescopes 
can be used to emulate six different two-telescope stereo systems: three combinations with 
intertelescope distances around 70 m (from the central telescope to the three corner telescopes), 
two combinations with about 90 m and 110 m (the sides of the imperfect square), and one 
combination with 140 m (the diagonal). Data from pairs of (triggered) telescopes are analyzed, 
ignoring the information provided by the other telescopes. Since only two telescopes are required 
to trigger the system, there is no trigger bias or other influence from those other telescopes. 

A slight difficulty arises since one compares combinations of different telescopes, rather than 
varying the distance between two given telescopes. While the four HEGRA system telescopes 
used here are identical in their construction, they differ somewhat in their age and hence the 
degree of mirror deterioration, in the quality of the alignment of the mirror tiles, and in the 
properties of the PMTs in the cameras, which show systematic variations between production 



batches. Different mirror reffectivities and PMT quantum efficiencies result in slightly different 
energy thresholds for given (identical) settings of electronics thresholds. The determination of 
image widths and hence the background rejection are sensitive to the quality of the mirror 
adjustment. These effects have to be determined from the data, and compensated. 

The analysis is based on the data set recorded in 1997 during the outburst of Mrk 501. Only 
data during the high-fiux periods between MJD 50595 and MJD 50613 were used. The data 
set was further restricted to small zenith angles, less than 20°, to approximate the situation for 
vertical showers. Mrk 501 was observed in the so-called wobble mode, with the source images 
0.5° away from the center of the cameras. An equivalent region imaged on the opposite side of 
the camera center was used as off-source region. Given the angular resolution of about 0.1°, the 
on-source and off-source regions are well separated. 



3 Analysis and results 



As a first step in the analysis, the detection rates were studied as a function of intertelescope 
distance. 

To equalize the energy thresholds of all telescopes, pairs of telescopes were used to reconstruct 
showers and events with cores located at equal distance from both telescopes were selected. By 
comparing the mean size of the images in the two telescopes, and in particular by comparing 
the mean signal amplitude in the second-highest pixel (which determines if a telescopes triggers 
or not), and can derive correction factors which can be used to equalize the response of the 
telescopes. Three of the four telescopes were found to be identical within a few %, one had 
a sensitivity which was lower by about 25%. In the analysis, pixel amplitudes are corrected 
correspondingly, and only camera images with two pixels above 20 photoelectrons are accepted. 
This software threshold is high enough to eliminate any influence of the hardware threshold 
(at about 8 to 10 photoelectrons, depending on the data set), even after the worst-case 25% 
correction is applied. 

The resulting detection rates of cosmic rays and of gamma-rays - after subtraction of the cosmic- 
ray background - are shown in Fig. 1. The rates measured for the three different telescope 
combinations around 70 m spacing agree within 5%, indicating the precision in the adjustment 
of telescope thresholds. Detection rates decrease with increasing distance; between 70 m and 
140 m, rates drop by about 1/3, with a very similar dependence for gamma rays as compared 
to cosmic rays. At first glance, this seems surprising, since compared to gamma-ray showers, 
the distribution of Cherenkov light in proton-induced showers - and hence the probability 
to trigger - is more strongly peaked near the shower axis, favoring a smaller separation of 
telescopes. On the other hand, however, the trigger probability for gamma-ray showers near 
threshold - where most of the detection rate originates - cuts off sharply around 120 m to 
130 m, whereas proton-induced showers occasionally trigger at larger distances [13]. The two 
effects compensate each other to a certain extent. In addition, about 33% of the cosmic-ray 
triggers are caused by primaries heavier than protons [14]. These particles interact higher in 
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Fig. 1. Detection rates of cosmic rays and of gamma rays as a function of the spacing between stereo- 
scopic pairs of telescopes, using a software trigger requirement of more than 20 photoelectrons in two 
pixels of both cameras. The normalization is arbitrary. Gamma-ray events were selected using very 
loose cuts on pointing relative to Mrk 501 (< 0.45°) and on the mean scaled width of the images (< 1.3), 
keeping essentially all gamma-rays. Cosmic-ray background in the on-source sample was subtracted 
on a statistical basis. The lines are drawn to guide the eye. 

the atmosphere and produce a wider light pool, again enhancing the rate for telescopes with 
wide spacing. In the later analysis, shape cuts completely eliminate such events. Indeed, if only 
cosmic-ray events with narrow (proton- or gamma-like) images are accepted, their rates fall off 
steeper with distance than the all-inclusive rate, or the gamma-ray rate. 



As a measure of the sensitivity for weak sources, the ratio S/y B of the gamma-ray rate to 
the square root of the cosmic-ray rate was used. The significance of background-dominated 
signals scales with this ratio. For optimum sensitivity, S/yB is optimized by tighter cuts on 
the pointing of showers, and on the image shapes. Fig. 2 shows the angular resolution provided 
by pairs of telescopes as a function of spacing. The three pairs at 70 m show differences at 
the level of 10%, indicating small differences in the quality of the mirror alignment and of 
the telescope alignment. Angular resolution improves slightly with increasing spacing; however, 
given the 10% systematic variations, this effect is of marginal significance. 

To determine the sensitivity for the different combinations of telescopes, cuts on pointing and on 
the mean scaled width were optimized for each combination, resulting in pointing cuts around 
0.1° and cuts on the mean scaled width around 1.05. In addition, a lower limit of 0.75 was 
imposed on the mean scaled width. Due to differences in the quality of the mirror alignment, 
the enhancement of significance due to such cuts differs slightly between telescopes. This can be 
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Fig. 2. Angular resolution provided by pairs of telescopes as a function of telescope spacing. A soft- 
ware trigger threshold of 20 photoelectrons is used. The angular resolution is determined by fitting 
a Gaussian to the difference between shower direction and direction to the source, projected on two 
orthogonal axes. The dashed line is drawn to guide the eye. 

studied by selecting events with telescopes A and B, and deriving the significance of the signal 
by cutting only on the event shape in A, or only in B. For identical telescopes, the results 
should be identical. Based on such studies, sensitivity corrections for the different telescope 
pairs were derived; these correction factors were always below 10%. Fig. 3(a) shows the resulting 
significance (defined as S/yB) of the Mrk 501 signal. The errors shown are statistical, and are 
dominated by the low statistics in the cosmic-ray background sample after cuts (even though 
the background region was chosen to be larger than the signal region). One finds that the 
significance is almost independent of telescope spacing, over the 70 m to 140 m range covered. 

Concerning systematic uncertainties, we note that the results for the three combinations at 
70 m agree reasonably well. Variations of parameters such as the software trigger threshold, or 
the exact values of the pointing cuts and angular cuts produce stable results, with systematic 
variations of at most 10%. 



The absolute significance of a point source detected with two stereoscopic lACTs, and its 
dependence on the intertelescope distance will, of course, depend on additional cuts which 
may be applied in the analysis. Additional requirements may, for example, concern the stereo 
angle, defined as the angle between the two views of the shower axis provided by the two 
telescopes, or, equivalently, the angle between the image axes in the two cameras. If the stereo 
angle is small - as it is the case for events with shower impact points on or near the line 
connecting the two telescopes, and for very distant showers - the two views coincide and the 
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Fig. 3. (a) Significance S/VB of the Mrk 501 signal for different pairs of telescopes, with optimized 
cuts on shower pointing and image shapes in the two telescopes. The significance is normalized to ~ 1 
for the 70 m spacing. The dashed line is drawn to guide the eye. (b) As (a), but with an additional 
cut of a stereo angle of at least 45° between the two views. The normalization if the same as for (a), 
so that the two data sets can be compared directly. A weaker 20° cut on the stereo angle, as used 
in some HEGRA analyses, has little effect on the significance, and yields the almost same results as 
shown in (a). 



spatial reconstruction of the shower axis is difficult. In some analyses, one will therefore add 
the requirement of a minimal stereo angle between the two views, in order to increase the 
reliability of the stereoscopic reconstruction. For a minimum stereo angle of 20° - as used in 
some HEGRA work - the resulting significance remains virtually unchanged compared to the 
data shown in Fig. 3(a). With a large minimum stereo angle of 45° (Fig. 3(b)), the distance 
dependence becomes more pronounced; for small telescope distances, significance is reduced, 
while for large distances it remains or is even slightly enhanced. The explanation is simple. For 
a telescope spacing d and a stereo angle greater than 45° , shower cores are accepted at most up 
to a distance of ~ 1.2d from the center of the telescope pair. As long as this distance is below 
the radius of the light pool, the effective detection area is reduced. One concludes that such a 
cut might be beneficial to reduce systematic effects in spectral measurements, for example, but 
should not be applied in searches for gamma-ray sources. 

Another quantity which infiuences the distance dependence is the field of view of the cameras. 
Cameras with small field of view will perform worse for large distances. For example, rejecting 
events with image centroids beyond 1° from the camera center reduces the sensitivity for the 
140 m point by a factor of about 2. 



4 Summary 



The performance of two-telescope stereoscopic lACT systems was studied experimentally using 
the HEGRA telescopes, as a function of telescope spacing in the range between 70 m and 
140 m. While detection rates decrease with increasing spacing, the significance of source signals 
is almost independent of distance, with a slight improvement for distances beyond 100 m. These 
results confirm Monte Carlo simulations (see, for example, [4]) which generally show that the 
exact choice of telescope spacing is not a very critical parameter in the design of lACT systems 
(provided that the field of view of the cameras is sufficiently large not to limit the range of 
shower impact parameters). An analysis requirement of a minimum angle between views favors 
larger distances. 

The results shown apply, strictly speaking, only to two-telescope systems. In lACT arrays 
with a large number of telescopes, one may place telescopes at maximum spacing in order 
to maximize the effective area of the array, under the condition that most individual showers 
are observed by two or maybe three or four telescopes. Alternatively one may choose to place 
telescopes close to each other, such that an individual shower is observed simultaneously by 
almost all telescopes, improving the quality of the reconstruction and lowering the threshold, 
at the expense of detection area at energies well above threshold. In the first case, one would 
probably use a spacing between adjacent telescopes between 100 m and 150 m; in the second 
case, one would limit the maximum spacing between any pair of telescopes to this distance. 
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